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Spectral densities are computed in unprecedented detail for quantum antiferromagnetic spin 1/2 
two-leg ladders. These results were obtained due to a major methodical advance achieved by opti- 
mally chosen unitary transformations. The approach is based on dressed integer excitations. Con- 
siderable weight is found at high energies in the two-particle sector. Precursors of fractional spinon 
physics occur supporting the conclusion that there is no necessity to resort to fractional excitations 
in order to describe features at higher energies. 
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Low-dimensional quantum antiferromagnets are of fun- 
damental and enduring interest. One reason is that high- 
temperature superconductivity depends crucially on the 
interplay of charge carriers with the magnetic excitations 
of a two-dimensional (2D) quantum antifcrromagnet. A 
vivid debate concerns the nature of the elementary exci- 
tations of such an antifcrromagnet. In terms of integer 
spin waves (S = 1), no quantitative description is avail- 
able of the spectral densities at higher energies, where 
a significant part of the spectral weight is locatedEm. 
Therefore, it has been suggested that fractional {S = 1/2 
spinons) excitations play an important role in 2D0 sup- 
porting the view that spinons and spin-charge separation 
are the basis of high-T c superconductivityQ'B. Spin lad- 
ders are similar to 2D planes in the range of higher en- 
ergies. The importance of fractional excitations for the 
description of high energy excitations in 2D would hence 
imply their importance in spin ladders. Here we present a 
description of the spectral densities in unprecedented de- 
tail in terms of integer excitations. Our results show that 
the essential point is not the fractionality of the elemen- 
tary excitations ("particles") but the proper description 
of multi-particle excitations. 

Spectral densities provide information on the density of 
elementary excitations, on their interaction and on how 
the particular excitation operator couples to them. For 
instance the dynamic structure factor as measured by 
inelastic neutron scattering couples to excitations with 
total spin S — 1. Integer spin excitations induce a dom- 
inant so-called quasi-particle peak in the dynamic struc- 
ture factora. If the integer S = 1 excitation decays into 
two fractional S = 1/2 spinons, which move indepen- 
dently at large distances, the quasi-particle peak vanishes 
and an important continuum appears at higher energies. 
The generic example is found in spin_chains for which 
spinons are the elementary excitations!]. 

We address 2-leg spin ladders, i.e. two chains with 
intra-chain coupling Jii coupled by an exchange coupling 
J± 

H = / ][J\\(8i t jSi t i+i + S2,iS2,i+i) + ■/±Si ] tS2,i] , (1) 

i 

where i denotes the rung and 1,2 the leg. As far as 



the excitations at higher energies are concerned, spin 
ladders constitute an intermediate system between ID 
chains and 2D square lattices. Additionally, they occur 
in a multitude of substancesa comprising also supercon- 
ducting systems!! 

What is the most effective way to describe the exci- 
tations? Two approaches are possible: (i) If the rung 
coupling J± dominates in ([!]), i.e. Hq :— H \ 7 _ Q the ex- 
citations are local triplets on each rung (rung-triplets). 
When J|| is switched on they start to hop from rung to 
rung. On increasing Jii the excitations become more and 
more extended; they are rung-triplets dressed by a mag- 
netically polarized environment, (ii) If the rung coupling 
Jl is weak it can be treated as perturbation which acts 
on the elementary spinons in both chains. The rung cou- 
pling binds two spinons on either chain thereby bind- 
ing two fractional S = 1/2 spinons to an integer 5=1 
triplet. The size of this bound object tends to infinity for 
Jl -> 0. 

For the value J± = Ju needed to compare with the 
square lattice it is a priori unclear which picture is supe- 
rior. Note that we are not aiming at the question which 
are the elementary excitations in the strict sense of the 
word since this question concerns only the lowest ener- 
gies. But we ask which excitations describe the physics 
at all energies best. We show that the approach (i) in 
terms of integral dressed excitations works excellently. It 
provides for the first time a quantitative description of 
the spectral densities on all energy scales. 

Previous investigations in gapped ID systems were also 
concerned with the nature of the excitations. For in- 
stance, Sushkov and Kotov compared quantum antifer- 
romagnetism to quantum chromodynamics in the sense 
that the 5=1/2 spins correspond to quarks and integer 
triplets to (vector) mesonsEl Zheng et al. described the 
deconfinement of fractional spinons in dimerized, gapped 
spin chains on decreasing dimerization in terms of integer 
triplet excitationaJ. 

We construct a mapping by a continuous unitary trans- 
formation (CUT)E3 of the Hamiltonian in (|I|) to an 
effective H c q conserving the number of rung-triplets: 
[Ho,H c ff] = 0. The ground state of H e g is the rung- 
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triplet vacuumEl. In this way a tremendously compli- 
cated many-body problem is reduced to a tractable few- 
body problem! The mapping needs an auxiliary vari- 
able I running between and oo. The Hamiltonian 
H(£) is transformed from I = (H(0) = H) to I = oo 
(iT(oo) = H eS ) by 

dH/dt=W),H{i)\ (2) 



Vij(£) = sign(Q M - Q S J )H,j{() , 



(3) 



where r\ defines the mapping. The matrix elements r\i 



parity with respect to reflection (RSc) about the cen- 
tcrlinc of the ladder implies an odd (even) number of 
rung-triplets. In Fig. 2 the distribution of spectral weight 
among the sectors of different number of rung-triplets 
is depicted. For local (i.e. fc-integrated) excitations the 
relative intensities P n — I n /Itot are plotted where the 
intensity I n stands for the w-integrated spectral densi- 
ties I(u>) in the sector of n rung-triplets. Recall that 
by our transformation the number of rung-triplets has 
become a good quantum number. The total intensities 



and Hi j are given in an eigen-basis of the number of ^tot = / + I(u)du = Y^Li In are accessible through the 



rung-triplets Q = Hq, Eqs. (|2j|3|) constitute a general, 
versatile prescription for any many-body problem to ob- 
tain an effective model in terms of elementary excitations, 
the number of which is counted by Q. The choice (^) 
eliminates all parts of H changing the number of triplets 
while retaining a certain simplicity in H(£) for interme- 
diate values of «L3o, namely that the number of rung- 
triplets is changed at most by ±2. For details on the form 
of H e g the reader is referred to Ref. [l3| To determine 
response functions the physical observable under study 
must be subject to the same unitary transformation (j2j) 
as the Hamiltonian. Hence the CUTs are based on a 
very clear-cut concept. Indeed, this concept renderecLthe 
computation of bound states in higher order possibleEjO. 
Conventional cluster techniquesO allow equally to deter- 
mine the energies of bound states very accurately if they 



are supplemented by similarity t 
tuting the unitary transformat 
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approach with hard-core bosonsr tr ! describes the excita- 
tions equally in terms of triplets with qualitatively similar 
results deviating, however, quantitatively for x 0.5. 

We were able to keep terms in H(£) up to order 14 
in x := J\\/J± for the kinetic energy (triplet hopping) 
and up to order 13 for the triplet-triplet interaction. Ob- 
servables were transformed up to order 10 allowing for 
the first time to determine true multi-particle continua, 
i.e. spectral densities, not only spectral weights, from a 
non-simulation approach. Our perturbative approach is 
supplemented by standard extrapolations and optimiza- 
tions of the resulting series. It turned out that it is pos- 
sible to determine the spectral densities within about 5% 
accuracy for x « 1. For lower values of x, the accuracy 
improves rapidly so that the results at x J$ 0.6 can be con- 
sidered exact. The results are not plagued by any finite 
size nor by any finite resolution effects thereby providing 
for the first time predictions ip-|great detail. 

For experimental relevanceEH and for comparison to 
the square lattice we present results for Ju = J±. For 
orientation Fig. 1 depicts the dispersion of the gapped 
elementary tripletsE3~Ej, the resulting lower and upper 
edge of the 2-triplet continuum and_thc. bound, states in 
the S = and the S = 1 charmeKBEl^'B The 
spectral densities are computed for the resolvent 

J( W ) = -tt"^ (0 \R(u + E - H)~ l R\ 0) rctardcd (4) 

for various operators R connecting the ground state to 
excited states of different spin and parity. Odd (even) 



sum rule 7 tot = (0|i? 2 |0) - (O^O) 2 . 

For x = 0, the excitation of a local rung-triplet ex- 
hausts the entire spectral weight for S = 1, i.e. I\(x = 
0) = 1. On increasing x the triplet becomes dressed 
and the quasi-particle weight decreases. The sum of all 
relative intensities must be unity. Up to x « 1 our re- 
sult (dashed lines) deviates from the sum rule less than 
3% giving evidence for the reliability of the extrapola- 
tions. For larger x the weight is overestimated increas- 
ingly. For x — 1 the 1\ and the 1% contributions exhaust 
74%+19%=93% for S = 1 and 77% for S = of the total 
weight. The neglect of contributions from more triplets 
is hence well justified. A description in terms of 1, 2 or 3 
integer excitations works perfectly without need to resort 
to spinons. 

In Figs. 3 and 4 the fc-resolved spectra are shown for 
n = 1 (Fig. 3a) and n = 2 (other panels). In both 
spin channels the 2-triplet spectral densities display im- 
portant and fine-structured continua bounded by square- 
root edges. Generally the interaction shifts weight down 
to lower energies. For momenta not to far away from 
7T, bound states (grey lines in Figs. 3 and 4) leave the 
continua. These bound states carry a large fraction of 
2-triplet weight. The differences between Figs. 4a and 
b result from the different parity of the two excitation 
operators at k — n with respect to reflection about the 
rung-direction (RSa), see also Ref. [H| 

A very intriguing feature is the line of square root sin- 
gularities seen within the continuum in Figs. 3b (dark 
grey line) and 4b. It increases from k ~ 7r/3 to k — n be- 
ing particularly clearly visible above k = ir/2. Below this 
line the larger part of the spectral weight is found. The 
momentum dependence of these mid-band square root 
singularities is strongly reminiscent of the upper limit of 
the 2-spinon continuum in single chainaZI, to which it will 
evolve for J± — > 0. The striking fact that a 2-triplet 
description is capable to yield 2-spinon features leads us 
to the conjecture that the important point is not the 
fractionality of the excitations but the correct descrip- 
tion of multi-particle excitations. Support for our con- 
jecture is provided by the computation p£ bound states of 
two spinonsE3'E3 in terms of two tripletsEj'EU in dimerized 
spin chains. 

In Fig. 4 for small values of A: a second upper peak 
occurs which is due to density-of-states effects resulting 
from the shallow dip that the 1-triplet dispersion (thin 
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solid line in Fig. 1) displays at k = 0. This local min- 
imum in the dispersion induces a van-Hove singularity 
which is smeared out by the interaction to the additional 
peak. It appears that the two parts of the dispersion 
from k — to k w OAtt and from k = ir to k ir/2 act 
as (almost) independent bands. With the help of this 
picture, the rather sharp resonance at about k = tt/2 in 
Fig. 4 is interpreted as bound state of the part of the 
dispersion for |fc| G [0, OAtt}. It is not an absolutely sta- 
ble state, i.e. not a <5-peak of zero width, because it may 
still decay into states with |fc| £ [0.47r,7r]. To be more 
specific, the dispersion for small values of k can be de- 
scribed approximately by w{k) k, Jj_(1.9 — 0.1 cos(fc/0.4)) 
(see Fig. 1) so that the two-particle kinetic energy 0(fc, q) 
at total momentum k and relative momentum q given by 
Vt{k,q) = uj(k/2 + q)+u>(k/2 — q) becomes dispersionless 
for k = 0.4vr: 0(0.471", q) w 3.8 J± for values of k/2 ± q in 
the range of the shallow dip. This implies that the rel- 
ative kinetic energy is negligible and interaction effects 
dominate ! Indeed, Fig. 4 shows that the resonance is 
shifted by attractive interaction from 3.8J^ to w 3.5Jj_. 

Our results predict many features in experiment. By 
inelastic neutron scattering the 2-triplet S = 1 response 
(Fig. 3b) is detectable if the momentum along the rungs 
can be tuned to zero so that only the symmetric com- 
bination (Sf { + S'| i )/2 leads to excitations. Then the 
2-triplet response (Fig. 3b) is separated clearly from the 
stronger 1-triplet response (Fig. 3a). To gain further in- 
sight in the 2-triplet interaction in the S = 1 channel we 
strongly suggest such experiments. 

Optical spectroscopy in terms of bimagnon-plus- 
phonon absorption detects a weighted superposition of 
the curves in Fig. 4 providing the first experimental 
verification q£,the existence of the bound state in the 
5 = septorEll. Raman spectroscopy measures in lead- 
ing ordeiO the curves at k = in Fig. 4e3. T-he theo- 
retical results agree very well with experimented though 
some deviations remain. These are are most probably 
due to the general presence .of paajuegligible cyclic ex- 
change terms in the cupratescJat^EJ. 

Our findings on the distribution of spectral weight 
show that large weight at high energies does not imply 
the necessity to resort to fractional excitations contrary 
to a frequently used line of argument, see e.g. Ref. |[ In- 
teger excitations allow to describe spectral densities of 
spin ladders in great detail on all energy scales. It turns 
out to be essential to describe multi-particle excitations 
adequately. We find also that this task is tractable since 
the sectors of low particle-number dominate clearly. 

We expect that a similar approach with integer excita- 
tions will also allow to describe undoped 2D cuprates 
quantitatively. Continuous unitary transformations 
proved to be a very powerful concept to deal with dressed 
excitations. The above findings call for a renewed discus- 
sion of the importance of fractional excitations in low- 
dimensional strongly correlated electron systems such as 
doped and undoped high-T c superconductors. 
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FIG. 1. Dispersions at Jii = J±. Elementary triplet 
(thin solid), lower and upper 2-triplet continuum edge (thin 
dashed), bound 2-triplet states with S = (thick solid) and 
S = 1 (thick dashed). 




FIG. 2. Relative intensities I r n of sub-sectors of different 
number n of triplets and their sum (dashed) for local excita- 
tions Ri. Left panel S = 1: Ri = Si^; right panel S = 0: 
flj = Si,iSi,i+i. 
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FIG. 3. fc-resolved spectral densities for operators indi- 
cated (S — 1); 5-peaks (grey lines) broadened for visualization 
by Ji/20 and inserted in front of the unbroadened continua. 
The bound state energies and the continuum edges are also 
shown (dashed lines), (a) 1-triplet peak ; (b) 2-triplet con- 
tinuum (multiplied by 4 to show details) and bound state. 
The dark grey line is a guide to the eye linking the mid-band 
square root singularities. 
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FIG. 4. Same as in Fig. 3 for S = 0; (a,b) 2-triplet contin- 
uum and bound state. In (b) continuum multiplied by 4 to 
show details. 
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